Creation and manipulation of anyons in a layered superconductor-2DEG system 
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We describe and analyze in detail our recent theoretical proposal for the realization and manipu- 
lation of anyons in a weakly interacting system consisting of a two-dimensional electron gas in the 
integer quantum Hall regime adjacent to a type-II superconducting film with an artificial array of 
pinning sites. The anyon is realized in response to a defect in the pinned vortex lattice and carries a 
charge ±e/2 and a statistical angle 7r/4. We establish this result, both analytically and numerically, 
in three complementary approaches; (i) a continuum model of two-dimensional electrons in the vor- 
tex lattice of the superconducting film; (ii) a minimal tight-binding lattice model that captures the 
essential features of the system; and (iii) an effective theory of the superconducting vortex lattice 
superposed on the integer quantum Hall state. We propose a novel method to measure the frac- 
tional charge directly in a bulk transport experiment and an all-electric setup for an "anyon shuttle" 
implementing the braiding operations. We briefly discuss conditions for fabricating the system in 
the lab and its potential applications in quantum information processing with non-Abelian anyons. 
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I. INTRODUCTION 

In two dimensions (2D) the wavefunction of indistin- 
guishable particles can have exotic exchange properties: 
upon the exchange of two particles known as "anyons" ^ 
it acquires a phase factor exp(i0) ^ ±1 with an arbitrary 
statistical angle 0. This is in contrast to the situation 
in higher dimensions where = or tt corresponding 
to bosons or fermions, respectively. The reason for this 
departure from the usual quantum statistics arises from 
the unique topological properties of 2D systems.— If the 
ground state of such identical particles is degenerate the 
exchange operations can be described by matrices in the 
degenerate subspace, which are generically non-Abelian. 
Such particles are referred to as non-Abelian anyons. It 
is also known that particles with fractional charge may 
arise, regardless of dimensionality, as the collective re- 
sponse of a many-body system to topologically nontrivial 
background fields 

The interest in anyons is not merely academic. The 
topological character of anyons means that local pertur- 
bations alter their properties exponentially weakly with 
a characteristic length of the order of the separation be- 
tween anyons. (For a single anyon this length scale is re- 
placed by the system size.) This "topological protection" 
has been argued to be useful for fault-tolerant quantum 
information processing)^ If the anyons have a rich enough 
non-Abelian structure, one may realize a universal quan- 
tum computer, in which computations are performed by 
braiding anyons.^ Otherwise, topological protection for a 
non-universal subset of operations may be obtained. In 
any case, the system can be considered, at least, as a 
topologically protected quantum memory. However, the 
insensitivity to local perturbations also poses a problem 
for manipulating anyons, which is necessary for braiding 
operations, and for the storage and retrieval of informa- 
tion from the quantum memory. 

How could anyons be realized in nature? Of course, the 
first condition is the restriction to two spatial dimensions. 



The fact that the world is three-dimensional immediately 
implies that anyons can only be realized as collective ex- 
citations of a many-body system whose constituent par- 
ticles are necessarily fermions or bosons. The emergence 
of such excitations constitutes the phenomenon of frac- 
tionalization in condensed-matter systems. 

Wilczebi^ proposed a simple model of an anyon as a 
bound state of a charge q and magnetic flux <I>, in which 
either the charge or the flux, or both, have a fractional 
value in units of the electron charge, e, and flux quan- 
tum $0 = hc/e, respectively. A closely related state is 
realized in a 2D electron gas (2DEG) in a perpendicu- 
lar magnetic field in the fractional quantum Hall (FQH) 
regime, described by the filling factor v. When 1/v is an 
odd integer the excitations carry a fractional charge ve 
(Ref. d) and have a statistical angle v-k (Ref. FQH 
physics is the canonical example of fractionlaization. The 
Coulomb interaction is understood to play an important 
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FIG. 1: (color online) The schematic diagram of the proposed 
device. A 2D electron gas is sandwiched between two type-II 
superconducting films. 
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role in stabilizing the FQH states, which may be consid- 
ered strongly correlated in the sense that they cannot be 
described by filling a set of single-particle states. 

One may naturally ask, then, if this is a necessary con- 
dition of fractionalization. It has been known for a long 
time that strong correlations are not necessary for frac- 
tionalization of charge in one and three dimensions ,^^1^ 
where exchange statistics are trivial. In two dimen- 
sions, fractional charge and statistics have recently been 
argued to arise in certain lattice modela^iii^'i^ii^ that 
preserve time-reversal symmetry and can be consid- 
ered weakly-interacting. Furthermore, we proposedi^ 
a weakly-interacting system to realize anyons that can 
be described by a Slater determinant of single-particle 
states. This proposal has the potential to be useful for 
manipulating anyons. Unlike some other proposals it also 
has a realistic chance to be fabricated in the lab. 

The system consists of a 2DEG in the integer quan- 
tum Hall state adjacent to a thin slab of a type-II su- 
perconductor. The idea is to employ the quantization 
of flux by the superconductor in units of 5 $01 i-e. half 
the natural quantum of flux of the 2DEG, as well as the 
quantization of Hall conductance in the quantum Hall 
state of the 2DEG. A sketch of a possible arrangement is 
shown in Fig. [T] We have chosen a symmetric placement 
of superconducting films around the 2DEG in order to 
minimize the spread and the in-plane component of the 
vortex magnetic field at the 2DEG. This is not necessary 
if similar conditions can be obtained in a different geom- 
etry. Both of these systems admit an effective single- 
particle description. The integer quantum Hall state 
can be understood as a full Landau level of essentially 
free electrons. The superconductor may be described 
as a condensate of Cooper pairs. It is then interest- 
ing that the proximity of these two weakly-interacting 
systems results in the fractionalization of charge and 
statistics. Anyons are formed in response to defects in a 
pinned Abrikosov vortex lattice in the superconducting 
film. Such defects may be detected and manipulated by 
scanning superconducting quantum interference devices 
(SQUID) and Hall bar probes or magnetic tips in mag- 
netic force microscopy that have been developed over the 
years i^^i^i^i^i^ Therefore, our anyons have the addi- 
tional potential of being manipulated despite their topo- 
logical protection. As we shall describe in the conclusion 
a similar setup with the FQH state at v = 5/2 should 
realize, and allow for the manipulation of, non-Abelian 
anyons. These anyons are robust as long as the spread 
of the vortex-lattice defect, controlled by the supercon- 
ducting penetration depth and the distance between the 
two layers, is small relative to the distance between the 
defects. 

In the present paper, in addition to providing the de- 
tails and some extensions of our previous analysiSfi^ we 
present new analytical and numerical work from alterna- 
tive starting points that shed light on different aspects 
of the system. Moreover, we present a novel and con- 
crete method to measure the fractional charge of anyons 



in a bulk transport experiment that utilizes the vortices 
in the superconducting film. We also describe an idea 
for the "anyon shuttle" , a system to manipulate anyons 
bound to the flux defects by purely electrical means. 

The rest of the paper is organized as follows. In Sec.HIl 
we describe the system and provide a general argument 
for our findings. In Sec. IIIII we discuss a continuum ap- 
proach to the electrodynamic response of the 2DEG to 
the vortex lattice and its defects, employing various sim- 
plifying assumptions and a combination of analytical and 
numerical methods. In Sec. IIVI we formulate and study 
a minimal lattice model of the system by exact diago- 
nalization numerically as well as in the continuum limit 
analytically. In Sec.|V]we study an effective theory of the 
system, which is expressed as a Chern-Simons-Maxwell 
theory that we use to study the interplay of the supercon- 
ducting and the topological orders in the system. Given 
the interesting physics predicted here and the potential 
application of the system in manipulating anyons, we be- 
lieve it is important to make an effort to fabricate it in 
the lab. Sec. IVII contains our proposal to measure the 
fractional charge and a description of the anyon shuttle. 
Finally, we conclude in Sec. IVIII by discussing the con- 
ditions for the experimental realization of the system. 
Details of some calculations as well as the derivation of 
some known results are given in Appendices EHCl in order 
to make the paper self-contained. 



II. THE SYSTEM 

The proposed system has two components: a layer 
hosting a 2DEG in an integer quantum Hall state, and 
a film of type-II superconductor. The purpose of the su- 
perconducting film is to quantize the magnetic fiux into 
an Abrikosov vortex lattice, where each vortex carries a 
fiux ^$01 which plays a central role in the physics we de- 
scribe. Cooper pair condensation in the superconducting 
film ensures that this quantization is extremely precise. 

As mentioned above, the anyon is to be realized in 
the 2DEG in response to a defect, i.e. an extra or miss- 
ing i$o flux in the vortex lattice. However, a naturally 
formed Abrikosov lattice is (almost) incompressible 1^ 
This means that if defects are to be created by mov- 
ing one of the vortices, the vortex lattice will rearrange 
itself in such a way as to compensate for the additional or 
missing fiux in the corresponding regions. Thus, in order 
to allow for the creation of such defects, we propose to 
artificially imprint an array of pinning sites on the su- 
perconducting film. The pinning sites are regions where 
superconducting order is weakened and can be created 
in a variety of ways.^^ They attract and pin supercon- 
ducting vortices, thereby preventing the incompressible 
rearrangement of vortices in response to a defect. At 
the matching field Bm the number of vortices equals the 
number of pinning sites. An increase (decrease) in the 
field away from Bm will induce a corresponding number 
of interstitial (vacancy) defects in the pinned vortex lat- 
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tice with a surplus (deficit) flux of ^$o- The interstitials 
can then be manipulated by a magnetic tipjSS, 

What is the response of the 2DEG system to a flux 
defect in the superconducting vortex lattice? To answer 
this question we will make two working assumptions. The 
first one is that the 2DEG is indeed in the integer quan- 
tum Hall state. In particular, we assume that the spatial 
variation in the field does not destroy this state. Second, 
we assume that the energetics of the system are dom- 
inated by the superconducting film. Importantly, this 
means that the quantization of fiux by Cooper pairs re- 
mains valid and exact. We shall discuss the conditions 
under which these assumptions hold true in Sec. IVIII 
Consequently, the question is reduced to: What is the 
response of the integer quantum Hall state in 2DEG to 
a surplus/deficit flux of 5^0? 

The answer is found by a thought experiment in which 
we slowly turn on the extra flux in time t from zero 
to 77<i>0 7 ^ = ij- The total charge accumulated at the 
position of the defect is 



dt CP Je 



n dZ, 



(1) 



where C is a contour in the plane containing the defect, n 
is a unit vector normal to it, and dl an element of length 
of C. In this process an electric field is induced in the 
plane by Faraday's law, E • dl = —d^/cdt. An electric 
field in the quantum Hall state with filling factor v results 
in a transverse electric current je = cthE x z, where z is 
the normal to the plane and an = ve^ /h is the quantized 
Hall conductance. Altogether, we find from Eq. ([T|), 



5Q — rji/e. 



(2) 



We note that for this result to be valid, the size of C 
needs to be much larger than both the size of the flux 
defect, and the size of the smallest Landau orbit. The 
former is determined by the penetration depth, Al, in 
the superconducting film, which comes in through the 
use of Farady's law. We shall see its role more clearly 
in Sees. IIII CI and |Vl The latter is of the order of the 
magnetic length £b = yhcJeB and is implied by the 
physics of the quantum Hall state, as will be discussed in 

secnni 

Thus, at = 1, we have uncovered a bound state of 
charge q — "qe and flux $ = Ty^g . This is an almost literal 
realization of Wilczek's model of an anyon.^° However, 
based on an analogy with the FQH quasiparticles, we ex- 
pect a statistical angle 9 = rfir = 7r/4, as opposed to 
the 277^7r — 7^/2 that would follow from Aharonov-Bohm 
phases. That this is the correct result can be affirmed by 
using the fusion rule'^ : the statistical angle of a compos- 
ite of n anyons with statistical angle 9 is rP'O. Putting 
two bound states (q, $) together we must obtain an elec- 
tron, a fermion: 40 = tt mod 27r. This is consistent with 
9 = tt/A. In effect, the anyon carries the memory of its 
fermionic past. 

These general findings will be confirmed by our de- 
tailed studies in the rest of the paper. 



III. CONTINUUM MODEL 

As a first concrete model of our system, we will study 
a continuum model of noninteracting electrons in the 
2DEG layer in magnetic field at filling factor v = 1. In a 
uniform magnetic field the energy levels are organized in 
Landau levels with a degeneracy iV$ = $/$o per spin, 
the total number of flux quanta in the plane. The ulti- 
mate goal of this study is to find the levels in the peri- 
odic magnetic field generated by the Abrikosov lattice of 
vortices in the superconducting film and in the presence 
of vacancies and interstitial defects. Not surprisingly, the 
full solution can only be obtained numerically, but we will 
find analytical solutions for simpler cases. Especially, we 
will use the general solution for the ground states of a 
2DEG in an arbitrary magnetic field, due to Aharonov 
and Casher^ in the special case of "Pauli" electrons, i.e. 
free electrons with the gyromagnetic ratio g = 2. 



A. Spin-polarized electrons 

We do not know the general form of the wavefunctions 
for other values of < g 7^ 2, but for large enough g ^ 1, 
the Zeeman coupling can be replaced with the condition 
that all electrons are in a single spin state, aligned with 
the magnetic field. We can then find the single-particle 
spectrum in the simple case of a uniform background 
field with vacancy and interstitial defects modeled by a 
(5-function profile. The single-particle Hamiltonian is 



1 



pol 



2mp 



P--A 

c 



(3) 



where p and A arc the momentum and the vector poten- 
tial operators. 

We shall work in the dimensionless polar coordinate 
{r cos ip,r sin If) — x/V2^b. The magnetic field is 
B(x) = {Bq + 77<I'o(5(x)) z. The vector potential is given, 
in the symmetric gauge, by 



1 $0 



(4) 



where 77 = ±i is the fraction of the flux quantum carried 
by the defect. 

The Schrodinger equation = E'i' is solved in Ap- 
pendixIXl From Eq. (jASp . the states in the lowest Landau 
level (LLL) are given by 



cx \z\ 



(5) 



with the complex coordinate z = re^"^. When 77 = 0, 
the LLL contains Nf, — N.^, states. For a vacancy, 77 < 
0, the orbital quantum number fc = 0, and the angular 
momentum I = 0, 1, • • • . All the states are pushed away 
from the center, so we still have A^e = Af$ states in the 
LLL and a charge deficit at the center. For an interstitial, 
77 > 0, on the other hand, k = 0, and / = 1, 2, • • • . All 
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the states are pulled in toward the center, so we lose the 
innermost state with I — 0, whose energy is pushed up 
into the gap. Therefore, the LLL now contains = 
A'^ — 1 states. The many-body ground state of the filled 
LLL is given by 



We 



(6) 



i=l 



i>] 



For further use, we also note a different form of this 
wavefunction in the "string gauge" where the vector 
potential of the additional (5-function flux is given by 
,5A = {^/2lBr)-^Ti<^o5{Lp)if. This gauge can be obtained 
from the symmetric gauge in Eq. (jH) by a gauge trans- 
formation A ^ A - VA, with A = (?7$o/27r)(p. The 
single-particle wave functions transform as ^' ^ ^'s = 
\I'exp(-i|^A) vpe-''''^. The effect of this gauge trans- 
formation is to send |zi| Zi in Eq. ([6|). 



N 



(7) 



i>j 



1. Fractional charge 

Using results from the previous section we have for the 
charge density in the thermodynamic limit 



i 

= e^|vl/|'')(x)p 



X.) 



where, 



1=0 



\l-ri\ 



T{l+\l^7j\ 



(8) 
(9) 

(10) 



The charge displacement due to additional flux can be 
calculated as SQrj = Q-q — Qo, where Qri — J p?)(x)dx = 
e F(T)dr. Note that Fq = 1. 
For a vacancy, ?7 < 0, 



FrXr) = 1 - 



r(-^,r) 



(11) 



where T{a, t) — J°° s° ''ds is the incomplete gamma 
function. For an interstitial, 77 > 0, 



F,M) = l- 



r(i-,7) r(i + ,7)- 



(12) 



In both cases, Fjj{t) ^ 1 for t 1, that is the length 
scale for the extra charge at the center is ^ £b- Thus, 



[F^{T)-Fo{T)]dT^r^e, 



(13) 



where we have used the identity r(a, T)dT = r(a-|- 1). 
This confirms the value of the fractional charge we found 
by our general argument in Sec. [IT] and clarifies the role 
of the magnetic length scale. 



2. Fractional statistics 

In order to find the statistics of the fractionally charged 
defects we calculate the Berry's phase^^ defined as the 
phase accumulated by a state "^w, as a parameter w is 
taken around a contour C, 



7B,w -T OB, in- 



(14) 



Here we have allowed for the possibility that '^^ is not 
an entire function of w (so it may depend on the complex 
conjugate w, too). When w parametrizes the encircling 
of one particle around the other, this is equal to twice 
the statistical angle. 

For a particle of charge q encircling an area with 
magnetic flux $, the Berry's phase is the same as the 
Aharonov-Bohm phase, 9ab = 27r(g/e)(<I>/<l>o). This 
property can be used to confirm the value of the charge 
we obtained above. The ground state in the string gauge 
for an extra fiux 77$o at a position with complex co- 
ordinate w is given by Eq. ([7]) with the replacement 
z~'^ {zi — w)~^ . In a large system of size L this is 
justified as long as \w\ <^ L / £b- A stronger argument 
will be given in the special case discussed in Sec. IIIIBI 
Denoting the ground state when 77 = as ^0(^1), we have 

vl/^ = ^^:^\z,,w) cc [](z, - w)-^^o{zi)- (15) 



So, 



' * m — 



(16) 



and dw'^w = 0. When w is taken around a contour C = 
dS this gives 



Of 



-i, / d^z / (^ll^^^fcill^d. 

Jc w ~ z 



= / d^z 



Pyjz - w) 
w — z 



dw 



(17) 



where we have used the definition ([9]) for the charge den- 
sity. Writing p,,, = po + Sp^^ we observe that the second 
term's contribution vanishes as with the size R of 
contour C (in units of £5). For the first term's contribu- 
tion, we note that the contour integral evaluates to 27ri 
when z G S' and otherwise, thus 



27r77 



O R- 



(18) 



where we have used the relationship between the charge 
and the magnetic fiux in 5, <I>s, for the v = I integer 
quantum Hall state. Comparing with 9ab we find the 
charge carried by the defects to be SQ = -qe as before. 
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The form of the wavefunction in ((T5)) suggests that we 
may write the wavefunction with two defects located at 
w and V as 



(19) 



We win justify this form more strongly in the next sec- 
tion. With this choice we may now compute the Berry's 
phase again, when w is taken around a contour C — dS 
encircling v. The algebra is completely analogous to the 
previous case and, when R — \v\ 3> 1, we find. 



2771] 

e 



(20) 



where the charge rje of the defect at v must now also be 
taken into account in p^. This gives 0b = 2nri{^s /^o+v) 
with the additional phase 29, 



(21) 



being twice the statistical angle, as expected. 



B. Pauli electrons 



Electrons carry a magnetic moment 



eh 1 

m„c 2 



(22) 



where cr is a vector of Pauli matrices, and couple to the 
magnetic field by a Zeeman interaction term, — /x • B. 
For free electrons the gyromagnetic ratio g ^ 2 with a 
very high precision. In a solid this value is renormalized 
and can be much higher. When g = 2 there is a power- 
ful method due to Aharonov and Casher^^ which gives 
the ground states of the free Pauli electrons moving in a 
general profile of a perpendicular magnetic field. 

The Hamiltonian for electrons in a perpendicular mag- 
netic field is 



2m„ ( 



P--A - fizB. 
c 



(23) 



The ground states can be found for an arbitrary space- 
dependent applied magnetic field i?(x) as tpoiz) ~ 
f{z)e~'^, where f{z) is an entire function of the (dimen- 
sionless) complex coordinate z and V^i/) = 27ri?/$o- For 
the sake of completness the derivation of this result is 
outlined in Appendix IB] 

We will now add a flux defect with a general profile to 
the uniform background, Bq, given by (5-B(x) = r/$o'5{(x) 
where 6^ is a general function that vanishes over a char- 
acteristic length ^ (in units oi Ib) and with unit weight 
/ ^^(x)d^x — 1. Then, switching to the dimensionless 
complex coordinate representation, we may write (/)(z) = 
i|2;p -|- ri(p^{z) where 



V20C(Z) = 27TS^iz), 



(24) 



with the solution given, using the Laplacian Green's func- 
tion (27r)^-'^ ln(r), as 



(j)^{z) = / \n\z — v\5^{v)d'^v. 



(25) 



For illustration, we note that when we take the limit 
^ — > we must have S(^{z) — s- d{z) and 4>^{z) ln|z| 
which reproduces the results of the previous section. This 
provides further justification for Eqs. and p^ . 

Transforming to the string gauge, the many-body 
ground state with a defect at position w is then 



{Zi 



w 



e-'"^«(^--'"%o(^«), (26) 



where refers to the ground state in the absence of the 
defect {f] = 0) as before. Thus, 



i 
i 

Using the expression 



dw(f>i{z 



1 1 



2w — Zi 



dw(f>i{z 



u — {z — w) 



2w — Zi 



,(27) 
.(28) 

(29) 



Similarly, 8^4'^ is found hy w w. Hence, follow- 
ing similar manipulations as in the previous section, the 
Berry's phase for transporting the defect around a con- 
tour C — dS is found to be 



/ d'zd'u 

2e 



Pviz - w;C) [S^iu) + S{u)] 



w 



(z — u) 



dw. 
(30) 

Note that the second term in the square bracket is exactly 
half of the result found in the previous section for ^ = 0. 
Similarly, Ob,w is found hy w ^ w and with [6^{u) — S{u)] 
in the numerator. 

Writing = po -\- Spri we see again that the correc- 
tion from the second term in 9b,w is 0(max[^, l]/i?)^, 
whereas the first term contributes exactly as before since 
/ d'^uS^iu) = J d'^uS{u) = 1. The latter also means that 
the leading contribution to 0b,w vanishes. Thus, 



27r77 



5s 



O 



max[^, 1] 
R 



(31) 



confirming the value of fractional charge, rye, as well as 
illustrating the role of the size of the defect, 

The statistical angle is calculated similarly from 9b 
by taking a defect at w around another defect fixed at 
V. The many-body wavefunction in the string gauge is 
given by 



^w.v{zi) (X 



-'^'^'('■''^-'"^^^{zi). (32) 



FIG. 2: (a) (color online) Charge density in units of e/2Txl% 
in a system with A^e = 100 electrons without (dashed blue 
line) and with (solid black line) a defect (j) = — i) of the 
widened flux tube with ^ = 0.5^s. (b) The difference in 
density —5prj = po — Pv The inset shows the accumulated 
integrated charge 5Q{r) = — / 5p^dx in units of e. 
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FIG. 3: (color online) Magnetic field (in units of Oo/tt^s) 
of the pinned Abrikosov lattice, (a) without, and (b) with 
a vacancy defect. The plots are a partial view of the larger 
system with 64 pinning sites, Xl = 0.2 and = 0.01 in units 
of vortex separation, ^/^i^B■ 



The algebra is the same as before, but now contains 
the charge of the defect at v so the leading contribution, 
when i? — |w| ^ max[^, 1], is given by ^-K-q^sl^^) + 26', 



(33) 



As an example, let us consider a widened flux tube, 
with a size ^ in units of \/2^_b, given by 



<5^(z) 



1 



which correctly tends to 5(z) as ^ — > 0. Then, 



(/)e(z)-ln/FF+^ 



(34) 



(35) 



which can be obtained directly from Eq. (p4|) . Thus the 
many-body wavefunction in the string gauge for a flux 
tube at w is found to be 



cx [](z, - iu)-" [1 + 



e 



(36) 



The derivatives needed for the Berry's phase calculation 
are 

d^^w = [ln(|2;, - +^2) + ln(z, - u;)] , 

(37) 

and similarly for di^^w These are the same as Eqs. (P7)) 
and (|28p. Finally, the Berry's phase is 



i— 



d^z i) Aw 



w — w — z 



(38) 

where = z — /{w — z). Therefore, when f <C i? the 
size of C = dS, both terms in Eq. ([55)1 have the same 
contribution, 2ti\ Jg(Pzp,^ = 27rie$s/$o + 0(l/i?^) to 
the integrals and we obtain Ob = 2nrj^sl ^o- The Berry's 
phase resulting from transporting one flux tube around 
the other can be calculated similarly to be 27r?7(<i>5/$o + 
T]) as long as max[^, 1] <C R. 



We have calculated the charge density and the accu- 
mulated charge in the vicinity of the widened flux tube 
numerically. The results are summarized in Fig. [2] which 
once again conflrms that the accumulated charge is ex- 
actly ie and is drawn from the edge of the system. 



C. Pauli electrons in Abrikosov lattice 

We will now use this general method to study the more 
realistic situation of defects in a pinned Abrikosov lattice, 
which we take to be square. Vortices are then separated 
by so that the flux through a unit cell (containing 

a vortex) is |$o- The average magnetic field is = 
(<i>o/27r^B)(l — Pd)-, where pd is the density of defects. 
The total field is obtained from the London equation 



(39) 



where denotes the position of vortices, is the coher- 
ence length, Al is the penetration depth, and is the 
profile of a single vortex in the plane. In the pure London 
model this would be a delta function resulting from the 
phase singularity at the vortex center. A more realistic 
approach takes into account the finite vortex core size 
which leads to a broadening of the delta function. The 
form that is easy to implement in calculations and also 
gives good agreement with the experimentally observed 
field distribution is given by a Gaussian^^ 



A«o(x) 



1 



exp 



2^0^ 



(40) 



If we let ^0 ^ we recover the ^-function vortices studied 
in Sec. [ifel 

We decompose the field B = Bq + SB + SB^^ where SB 
is the periodic modulation by the Abrikosov lattice and 
SBd is the field produced by defects. We choose to work 
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FIG. 4: (color online) The electronic density (in units of e) in 
the Abrikosov lattice (a, c) without, and (b, d) with defects. 
Panels (c) and (d) are zooms on (a) and (b) where a vacancy 
is introduced. The difference between (a) and (b) is shown in 
(e), and between (c) and (d) in (f). The system has 64 pinning 
sites, Xl ~ 0.2 and = 0.01 in units of vortex separation. 



in the Landau gauge where (j) = 2/^/2^^ from Eq. (|B6|) . 
We have returned to the dimensionful quantities for clar- 
ity. This corresponds to periodic boundary conditions in 
the X direction. Accordingly we choose the complete set 
of entire functions, fi{z) = exp[— 27riZ(a; + i?;)/Lj;], where 
I G iNl and Lx is the size of the system in the x direc- 
tion. The (un-normalized) single-particle ground-state 
wavefunctions arc then given by 



(41) 



where V^^^ = |^(5B and V^fJ^^ = ^^-S^. 

Using the periodicity of the system we have SB{x.) 



iGx 



SB. 



1 c2/^2 
" 2 ^O'-' 



G 



1 



(42) 



Here G is a reciprocal lattice vector. The G = is ex- 
cluded since the average extra flux J 6BdS = 0. Fig. [3] 
shows the magnetic field, from Eq. ((42|) . in the pinned 
Abrikosov lattice when ^ ^ ^b- This is the limit 
we are interested in: the first part ensuring that the su- 
perconductor is a strong type-II and the vortex core is 
small, and the second part ensuring that the vortices are 
well separated and defects are well localized. The result 
is a strongly modulated magnetic field. 

The Laplace equation (|B2p can be solved in the recip- 
rocal space, (5(/)(x) = ^ 



i'Ge 



iG X 



CG = 



SB, 



1 c2/^2 



G 



£|G2 1 



(43) 



We define a supercell of Nd cells of the original lattice 
that contains one defect; thus, pd = ^/N^- This super- 
lattice of defects has a corresponding reciprocal superlat- 
tice, whose vectors we denote by G^- Now, 



6d(x) ^ Pd ^ c 



OG^e 



iGd-x 



(44) 



Eigenfimctions ^'J(x) in Eq. (|4ip are linearly indepen- 
dent but in general do not form an orthonormal set. This 
complicates computation of observables. In the next step 
we thus map the states in (PT|) to an orthonormal basis 
by diagonalizing the overlap matrix 



(*;i*L> 



(45) 



Denoting the orthogonal matrix that does so by U, 
{U^AU)ik = aiSik with eigenvalues a;, the orthonormal 
basis is given by 



1 



k 



/a, 



(46) 



Then, we may easily obtain the charge density in the 
many-body ground state. 



P(x) 



I 



(47) 



We have employed Fast Fourier Transform methods— 
to obtain the charge density in the system when defects 
are introduced in the vortex lattice. In Fig. 2] we show 
our results for the charge density. In the regular pinned 
vortex lattice, the electrons are mostly uniformly spread 
on the plane with regular peaks bound to the vortices. 
At = 1, the average density is half an electron per 
vortex. A vacancy depletes the charge from around the 
vacant site, with an integrated charge deficiency which is 
exactly ^e. This result allows us to conclude that a vor- 
tex defect binds exactly quantized fractional charge even 
when the periodic structure of the underlying magnetic 
field is taken into consideration. 
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(a) (b) 




FIG. 5: (color online) (a) The Peierls factors, e'^^'J , in the 
Landau gauge: (green) dashed links are —1, (black) solid links 
+1, and (blue) diagonal links +i in the direction of the arrow, 
(b) The string gauge for a defect (black disk) : every link that 
intersects the string (dashed wiggly line) acquires an extra 
-1. 



IV. LATTICE MODEL 

In this section we consider a lattice model to study the 
response of the 2DEG layer to the arrangement of vor- 
tices in the superconducting film. The continuum model 
in the previous section and the lattice model in this sec- 
tion can be thought of as the same system for different 
realizations of the gyromagnetic ratio. For large values of 
g we can think of electrons as tightly bound to vortices. 
In this limit we can model the 2DEG with a tight-binding 
Hamiltonian for spin-polarized fermions on a square lat- 
tice where each plaquette is threaded with half a flux 
quantum, 5<&o- 

The Hamiltonian is, 

Hiatt = ^^J-e"'" clc^. + e,clq, (48) 

where a annihilates an electron at site Xi, a is an on- 
site chemical potential, and e^^'^ are (gauge-dependent) 
Peierls factors incorporating the magnetic flux: 

= -j- / A • dx, (49) 

with A the vector potential. The simplest model that 
produces the integer quantum Hall state is one with near- 
est and next-nearest neighbor hopping, t and t' , respec- 
tively. The next-nearest neighbor hopping is needed to 
break the time-reversal symmetry. When t' — the 
smallest flux an electron sees is through a plaque- 
tte, which is changed under time-reversal operation by 
$0- Therefore, the Hamiltonian is time-reversal invari- 
ant when t' = and the Hall conductance vanishes. 

For a flux i$o uniformly distributed over a plaquette 
the vector potential in the Landau gauge is given by A = 
|sgn(i?)<I>o(0, x), where B is the field perpendicular to 



the plane. The Peierls factors are then given by 

Xi,i+x = 0, Xi,i+y = T^ix, (50) 
Xi,i+x+y = Xi+x,i+y = Tr{ix + ^sgn(B)). (51) 

These Peierls factors are pictured in Fig. [5]Ja). We may 
then choose a two-site unit cell + x) and form the 
spinor field ipi — {ci, Q+s)"^. For a uniform onsite poten- 
tial Ci = e the Hamiltonian can be written compactly in 
the reciprocal space, with the reduced Brillouin zone {k : 

\K\ < TT/2,\ky\ < tt}, as Hutt = Ek'Al(e-2ii/k)?/'k, 
where 

i?k = -h{k)c7^ + ReD{k)c7^ + lmD{k)c7y. (52) 

Here crx,y,z are the Pauli matrices, and 

ft-(k) = cosky, (53) 
£'(k) — e'''^ (cos fca; — imsin fcj, sin /cj,) , 

where we have introduced the "mass" m — 2sgn{B)t' /t. 
The energy spectrum is then given by 



£'(k) = e ± 2t Y cos^ + cos^ ky + m? svc? k^ svc? ky, 

(54) 

which is symmetric around e and has a gap 8t' at the two 
independent "nodes" = (^,±-|). The symmetry of 
the spectrum is a general property of Eq. (|52p. since 

ayH^cTy = -iJk- (55) 

Therefore, CTytp'^ is an eigenstate of energy e — E, if V'b 
is one with energy e + E. The spectrum and the corre- 
sponding density of states are shown in Fig. [SJa-d) . 

A defect is introduced by an additional ?7$o flux 
through one of the plaquettes. This will alter the Peierls 
factors by Sxij, so that J^y'^Xy = ^ mod 27r around 
a closed loop containing the defect, and zero otherwise. 
We choose to work in the "string gauge" specified by a 
string originating from the defect and ending at a bound- 
ary: 5x(ij) = 27r77 if the string cuts the bond Xj — x^, and 
zero otherwise. This is shown in Fig.[51[b). Two different 
choices of the string are related by a gauge transforma- 
tion. It is important to note that our lattice formulation 
does not distinguish between an interstitial {rj = +i) 
and a vacancy {r] = — ^), since the difference in fiux is 
a full flux quantum through the smallest loop of the lat- 
tice (i.e. half a plaquette). As already mentioned above 
lattice electrons cannot distinguish such fluxes. 



A. Hall conductance 

At half filling and for t' Q the lattice model exhibits 
precisely quantized Hall conductance aH = ±(e^//i). To 
see this we employ the TKNN formula^* which gives 
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Oh = /C(e^//i) in terms of the integer topological invari- 
ant 



with 



JC = ^ (f dk - Ak, 

JdBZ 



-4k = ^ (V'skliVkV'sk)- 



(56) 



(57) 



Here ipsk are eigenstates of with band index s and 
the sum is over occupied bands. 

The integral in Eq. ([55]) represents the total Berry 
flux through the Brillouin zone. The largest contribu- 
tion comes from the vicinity of the nodal points and it 
is easiest to evaluate JC by first linearizing near the 
nodes and then computing their contributions separately. 
As shown below in Sec. lIVDI the linearized Hamiltonian 
near the node at has the generic Dirac form 



hp = PxCTx + PyCTy + m<7z 



(58) 



with the spectrum Ep = i-v/p^ + m'^ and p is the mo- 
mentum relative to the nodal point. A quick analysis of 
Eqs. ([55]) reveals that the mass terms at the two nodes 
are equal. A straightforward calculation summarized in 
Appendix [C] yields 



IC^ = isgn(m). 
Adding the two contributions, we thus find 



(7H = sgn(B) 



h ' 



(59) 



(60) 



as expected. This type of calculation is valid for any 
system with Dirac nodes. In graphene, for example, the 
mass terms created by a simple charge density wave have 
opposite sign and the Hall conductance vanishes, consis- 
tent with the fact that the system is time-reversal invari- 
ant in zero applied field. 




distance 



20 40 60 80 

separation 



FIG. 6: (color online) The spectrum of the lattice model (a) 
without a gap, t' — 0, and (b) with a gap t' /t — 0.1. The 
density of states in the (c) gapless, and (d) gapped system. 
In all cases e = is assumed. Charge density with (e) a 
single defect, and (f) two defects in a 36 x 36 system and 
t' /t — 0.3. (g) The integrated charge in a 40 x 40 system with 
open boundary conditions, (h) The energy splitting of zero 
modes of two defects vs. their separation in a 164 x 10 system 
with periodic boundary conditions. 



B. Exact diagonalization 

We have performed exact diagonalizations the lattice 
Hamiltonian ([^5]) on lattice sizes up to 50 x 50 in vari- 
ous settings. The system with a single defect and open 
boundary conditions supports a zero-energy bound state. 
Since the spectrum remains symmetric, standard argu- 
ments^^ then show that the charge bound to the defect 
must be ±56 at half-filling depending on whether the 
zero mode is filled or empty. As seen in Fig. [61(g) the 
charge bound to a defect indeed integrates to with a 
numerical precision within machine accuracy. The charge 
density profile in the ground state is shown in Fig.[6{e,f). 
The extra charge is localized around the defect center 
within a lengthscale ~ m~^. For two defects, the two 
zero modes are slightly split due to tunneling between 



them with the energy splitting that decays exponentially 
with the defect separation. This splitting is plotted in 
Fig. IDIh). The effects of disorder in the on-site potential 
ei were discussed in Ref. Specifically, we found that 
precisely quantized fractional charge persists even in the 
presence of disorder, as long as it is weak in comparison 
to the excitation gap. 

C. Lattice Berry's phase 

Since we have the complete spectrum of the lattice 
Hamiltonian, we can also calculate the Berry's phase ac- 
cumulated by the ground state of the system with two 
defects as one is taken around the other. 
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Consider the ground state \w) with two defects, one 
placed at a fixed position and the other at w. As we take 
the second defect around the first one through the posi- 
tions wi, W2, ■ ■ ■ , wn = wi, we may calculate the accu- 
mulated phase at step n though a generalized Bargamann 
invarianti^i^'^ 

Ob.u = arg(wi|u;2) (W2IW3) • ■ ■ {wn\wi) . 

(61) 

In each step the phase of the overlap changes by an in- 
cremental amount; the phase of the product of overlaps 
is then equal to the sum of all such incremental changes. 
Of course \w) is defined up to an arbitrary phase. The 
product in is independent of this phase for all states. 
Especially, the last overlap {wn\wi) is included to make 
the product independent of the arbitrary phase of the 
initial and final states as well. 

There is also a local gauge freedom in the Hamil- 
tonian (|48p and thus in \w). Therefore, the overlaps 
{wr\ws) in Eq. ((6T|) are gauge-dependent. This can be 
avoided by choosing to work with a gauge-invariant over- 
lap. 



{Wr\Ws) 



{wr\e-'^°'"\ws), 



(62) 



where xo,rsi^i) — ^ J^^ ^o(^,Xi)dt, t is time, and Aq is 
the electromagnetic scalar potential in the chosen gauge. 
Basically, Xo.rs(xi) is the temporal Peierls factor at site 

The significance of xo can be understood by think- 
ing about the dimensionless gauge- invariant flux $(C) = 
f^A-dx for a temporal loop Cij^rs = (ir, Xi) (t^, x^) 
(ts,Xj) {tr,Xj) — > {tr,Xi). On the one hand. 



<f(c.„..) 



$0 



2tt 

On the other hand. 



[XO,rs(Xj) - XO,rs(Xj) -|- Xijits) " XtjiU)] ■ 

(63) 



E • dx, 



(64) 



where E — — VAq — -^^tA is the in-plane electric field. 
This is simply the total integrated electromotive force 
along the line connecting to Xj. For each flux quan- 
tum crossing x^ — x^ in the time interval tg — tr, we see 
from Eq. ([M)) that ^{Cij^rs) changes by ±27r, coming 
from the electric field induced by the motion of the flux. 
Therefore, xo in Eq. ensures that this induced field 
is correctly taken into account. 

We shall work in the temporal gauge where Aq = 0, 
where the gauge-invariant overlap is the same as the reg- 
ular overlap. For a moving defect in the string gauge, 
this is satisfied when we extend the Dirac string exactly 
along the path of the defect. This can be seen by fo- 
cusing on a primitive temporal plaquette, p = C(y)^(rs)- 
If the defect crosses the spatial side of p, $(p) = ±7r. 
Since the Dirac string trails along the path of the defect, 
it also crosses the same spatial side of p, contributing 




FIG. 7: (color online) The temporal string gauge. The verti- 
cal direction represents the time and style/color is as in Fig.[S] 
As the defect is moved from a plaquette to the neighboring 
one, it threads the shaded temporal plaquettes, thus inducing 
an electric field. In the next time step the string is extended 
along the defect's path, which accounts for the induced elec- 
tric field. The Peierls factors in the vertical direction remain 
equal to +1. 
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FIG. 8: The Berry's phase, 6b /t^, and even-odd effect. In (a, 
b) there is a single defect in the system that is taken around 
a loop with (a) an even, or (b) an odd number of sites. In 
(c, d) the system contains two defects, of which one is taken 
around the other on a loop containing (c) an even, or (d) an 
odd number of sites. 



exactly ±7r to $(p) through the spatial Peierls factors 
in Eq. ([55)1 . Thus, ao^(,s)(xi) = ag (rs)(xj) — 0. If the 
defect does not cross p, nor does the Dirac string. This 
is depicted in Fig. [T] The same conclusion can be drawn 
in the continuum formulation. In the string gauge A is 
perpendicular to the Dirac string, therefore in order that 
Aq = A ■ v/c = 0, the velocity, v, of the defect must 
be tangent to the string, i.e. the Dirac string should be 
extended along the path of the defect. 

Fig. [5] summarizes our numerical calculation of the 
Berry's phase. When a single defect is taken around a 
loop C = dS, the net effect can be understood as the 
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Aharonov-Bohm phase of transporting the flux of the de- 
fect \r]\^o around the total charge inside the loop l-ZVg, 
where Ns is the number of sites enclosed by C. Hence, 
Ob = ^{Ns mod 2). This "even-odd effect" is seen in 
Figs.[8l^a) and[8l^b). When a second defect is introduced 
inside the loop C, the charge is reduced by |, which is 
equivalent in effect to reducing the number of sites by 1. 
Thus the "even-odd effect" must switch, as depicted in 
Figs. [DJc) andlUd). The statistical angle can be inferred 
from the difference in the Berry's phase with and without 
the second defect, giving 29 — tt/2 as expected. 

D. Low-energy limit 

Fractionl charge and statistics can also be found an- 
alytically by studying the Hamiltonian ([48]) at low en- 
ergies. This will give us a continuum description of the 
lattice model valid for the long-distance, collective be- 
havior of the system. 

At half-filling the low-energy spectrum is dominated 
by quasiparticles around the two nodes K^. Expand- 
ing around the nodes, k = + p, the low-energy 
Hamiltonian is found to be Tiiatt ^ J2p ^p^p^p^ with 
Tfp = (.i^K++p,i'K-+p), and 

Hp = aipx + a2Py + imaia2, (65) 

where a — (ai, 02) — — (I^CTj,, CTz^CTz). Since these ma- 
trices are diagonal in the nodal index, the two nodes are 
decoupled in the low-energy approximation. Moreover, 
by a rotation S = exp(i7rl (g) crz/4) exp(i7rl (g) o-;r/4)(l © 
ax) we find Hp HpS = hp, where hp is given 

by Eq. dMl). 

The effect of additional flux is to shift the momenta 
p ^ p — f (5A, where SA is the vector potential of the 
additional field. For an additional flux rj^o in the sym- 
metric gauge, 6 A = (r/<I>o/27rr)<^. In real space, then, 

H = a ■ (p — '^-^^ + imaia2, (66) 

where now p — — iftV is the momentum operator. 

The problem is reduced to Dirac fermions interacting 
with a (vanishingly thin) solenoid of flux ij^q. It is known 
that the spectrum of this problem is not completely spec- 
ified without extra boundary conditions at the origini^ 
Mathematically, this is related to having more than one 
self-adjoint extension of the Hamiltonian These 
self-adjoint extensions can be parametrized by a pair 
of angular variables {9^^j,9~^^), one for each node. For 

tt/2 < 9f^- < 37r/2 there is a bound state for the ± node 
with energy \E\ < m at any value of rj. When 9^^- = n 
the bound state at 77 = ±i has zero energy. The bound 
state wavefunction in each case is square-integrable but 
diverges at the origin. We also note that without fixing 
^adj actually finds a continuum of bound states. This 
is due to the scale invariance of the Hamiltonian (pS)) . 
which means that by scaling the distance r Xr we can 



go from a bound state ^hir) of energy E to another one 
^'b(Ar) with energy E/X. This anomalous bound state 
spectrum collapses onto a single bound state energy by 
properly regularizing the problem at short distances. The 
self-adjoint parametrization provides such a regulariza- 
tion procedure. 

The lattice itself provides a natural short-scale regular- 
ization, and our numerical results above show that indeed 
we have a single zero mode for each value of jryj = ^. 
Standard arguments^ then lead to a fractional charge 
SQ = rje. The sign is determined by noting'^'^ that the 
bound state energy tends to zero from negative (positive) 
values when 77 is tuned from to (— |) and therefore 
is to be included (excluded) in the ground state at half- 
filling. 

To find the statistics we must go beyond the spec- 
tral properties of the low-energy theory to an effective 
field theory description. This can be done in a path- 
integral formulation where ip takes on the role of a dy- 
namical field coupled to the electromagnetic gauge field 
SA — {SAo,SA) of the defects. The mass term \ma1a2 
breaks the time-reversal and parity symmetries in the 
Dirac Lagrangian. Integrating out fermionic field T at 
the one-loop order straightforwardly we obtain the effec- 
tive action S[5A\ ~ f dtdx£cs + ■ ■ ■ , where 

Ccs = '^^e'^^'SA^d.SA,, (67) 

is the Chern-Simons (CS) Lagrangian describing the 
topological sector of the effective theory and the dots 
indicate the remaining (non-topological part) of the the- 
ory. Of course, it is well-known that a CS term describes 
the topological aspects of QHE. The only difference here 
with the usual CS theory of QHE is that since the gauge 
field Sd describes the dynamics of defects in the (pinned) 
supercondcuting vortex lattice, it carries a magnetic flux 
quantized in units of ^$0 instead of ^o- 

The charge and statsitics of a defect carrying flux ?7$o 
can now be calculated in the standard way'^- (see also 
Sec. IV]) from the CS theory, yielding the charge SQ = rje 
and the statistical angle 9 = 77^ tt, as obtained previously. 



V. EFFECTIVE THEORY 

In this section we formulate an effective theory to de- 
scribe the proposed system and the interaction between 
the 2DEG and the superconducting film. Using standard 
duality mappings, we show that the topological proper- 
ties of the system in the vortex lattice phase are described 
by a Chern-Simons-Maxwell theory. Using this effective 
theory we derive the fractional charge and statistics of 
its excitations. 

Our notation in this section is as follows: the space- 
time vector is denoted by x = (a;'^,x^,x^) with x = 
{x^^x^) the spatial and = ct the temporal parts. 
We use the relativistic notation in Minkowski space-time 
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u ■ V 



L^^v" and (9 x uY = e^^^dyUx, where 



g — diag(— 1, 1,1) is the Minkowski metric. 



Formulation 



The 2DEG at odd inverse filhng factor l/v\s described 
by the Chcrn-Simons effective action^ 



'S'cs — 



1 ^ e - ^ ' 

a • (9 X a) H a ■ {d y. A) 

Attv 2tt 



(68) 



We remind the reader that the topological gauge field a 
couples to the electromagnetic gauge field A in just the 
right way as to lock the electromagnetic field d x A to 
the electronic current je — SScs/SA ~ ^{d x a): 



SScs 

Sd 







je ~ i^e d X A. 



(69) 



In particular, the electron density = veB. 

We shall describe the superconducting film at zero tem- 
perature by the action S'pin + Sxy where S'pin represents 
the effects of the lattice of pinning sites and S'xy is the 
3DXY model, which we first write on the discretized Eu- 
clidean space-time. 



S'- 



XY 



£e^A). 



(70) 



Here, kq is the superconducting phase stiffness, is 
the superconducting phase, e* = 2e is the Cooper pair 
charge, A is the lattice difference operator and the sum 
is over the points of a space-time cubic lattice with the 
spacing £ explicitly included to make the transition to 
the continuum limit transparent. This action describes 
the phase ordering of the superconducting transition in 
a magnetic field with the London penetration depth 
Xl = 1 / •\/47rKoef . In the continuum limit it could be ob- 
tained from the Ginzburg-Landau theory by taking the 
amplitude of the superconducting order parameter to be 
constant. Using the Villain approximation and standard 
duality mappings'^^''^® we may rewrite the 3DXY action 
in Minkowski space-time continuum in the "current rep- 
resentation," 



S'xy — 



d^a 



-^{dx s)^ -e.idx A) 
2ko 



(71) 

where s is the dual vortex field and jy is the conserved 
vortex current, d ■ jy = 0. 

In the following we shall take the main effect of Spin 
to be the pinning of the background vortex lattice and to 
provide a one-body potential for the defects that pref- 
erentially positions them at the interstitial locations. 
Therefore we will drop this term from the dynamics with 
no change in topological properties. 



B. Phases of the system 

The 3DXY action ([7T|) has two phases: (1) a vortex- 
lattice phase, where (s^) — in the dual or (e"') 7^ 
in the original problem; (2) a vortex-liquid phase, where 
(s^) ^ in the dual or (e''') = 0. We wiU now derive 
the effective action of the combined system. 



S, 



cs 



5XY, 



(72) 



in each of these two phases. 

(1) In the vortex-lattice phase we may integrate out 
the dual vortex field s by noting that if we write the con- 
served current jv = -^d x (j) where is a hydrodynamic 
fiux variable, we have 



1 



-77^(5 X s)2 -f (0 e,A) -{Oxs) 



SxY = J d'x 
Since (s^) = we may integrate s out to find 



SxY 



2 



d^x(0 - e^A) 



2 

±1 



(73) 



(74) 



where the subscript _L indicates that only the transverse 
parts of the fields enter this term. At this point we sep- 
arate the static and dynamic parts of A and (/> coming 
from, respectively, the pinned background vortex lattice 
and the defects: A — Aq + SA and <j) = (po + dcj). Then 
we may shift the electromagnetic gauge field of defects 
SA 5A + 5(j)/ei, and then integrate out 5 A to find the 
Chern-Simons-Maxwell form 



Scff — 



-Jv 



1 

( 



{d X a) 



1 

2k 



{d X a) 



(75) 

where k — (27re*/e)^Ko and jy is the total vortex current 
of the pinned lattice and defects. We shall study the 
topological properties of this action shortly. 

(2) In the vortex-liquid phase, the dual field s acquires 
a mass and the Maxwell term is replaced with a mass 
term 



S 



XY 



d^a 



1 

2^0 



+ s-dx{$-e,A) (76) 



d X e^A) 



(77) 



after integrating s out. Shifting and integrating out 6 A 
as before we now find a mass for a: 



d^x 



—3v ■ a - a ■{dxa) + -a^ 



(78) 



where /i = (e/27re^)'^/io. Thus in the vortex liquid phase 
the gauge field a is gapped and its topological effects are 
absent. 
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C. Fractional charge and statistics 

The fractional charge of defects in the vortex lat- 
tice phase, i.e. the Chern-Simons-Maxwell effective the- 
ory (|75p . is found as follows. From Eq. ([75)1 we find 
SSes/Sa — gives 

3e~^dxje- \ri\iyejy = 0, (79) 

where — x a, jryj = e/e* = i, and ^ = 27ti'k^^ is 
the coherence length. This modifies Eq. ([69|l . reflecting 
the effect of the superconducting film. After a Fourier 
transform we find 



1 - 



jv + i^fc X jy - ^'^{k-jy)k 



(80) 



So, the excess charge density vanishes exponentially over 
a distance ^, which is basically the size of the defect. The 
total charge Qit) = J dxf{t, x) k = 0), which by 

Eq. (EOl) is 



d^x 



27r 



(a X ,/.)"', 



(81) 



and we have used the definition of the vortex current in 
terms of the flux variable (j). The integral here is noth- 
ing but the total vorticity (the winding number) of the 
lattice, a topological invariant. Introducing an intersti- 
tial (vacancy) at the origin changes the vorticity by -1-1 
(— 1). Thus, defining the sign of 77 to be the same as the 
vorticity, the charge bound to a defect is found to be 



SQ = rjiye. 



(82) 



The fractional statistics can be established by a Berry's 
phase calculation in the effective theory. To this end 
let's write the vortex current in Eq. ((75)) as jy — jo + 
Sji+Sj2 where jo is the vortex current of the background 
pinned lattice and 6j^{t,x.) = {dx'^ /dt)vsS{x — Xs{t)) is 
the vortex current for defects s = 1,2, with world-lines 
Xs = (ct,Xs) and vorticity Vg — ±1. We take a contour 
r2 for X2 (t) that encircles the stationary position xi = 
and ask what is the Berry's phase contribution from the 
cross terms of Sji and Sj2- 

By integrating out the Chern-Simons field we find 



In order to invert the matrix M in Eq. ([84)) . we need 
to fix the gauge for a. We do this by adding a term 
Mjf — > M^^+Qk^^ky. In general C, can also be fc-dependent 
producing a nonlocal gauge fixing term. The result is 



Af-i = ^ 

A"^ 1 - ^2/,2 
1 



Cfc2 1 - ,f2/j2 



(86) 



By current conservation, d-js — 0, the term k^k^ does 
not contribute to 9b, as is required by gauge invariance. 
After some tedious but ultimately standard calculation 
we find 9b — 9dyn + ^top, where the topological (dynam- 
ical) phase is contributed by the off-diagonal (diagonal) 
part of AI^^ in Eg. ([55)1 . The topological phase is twice 
the statistical angle, 6'top = 26*. 
The dynamical phase 



^dyn = vriiri2- I dtKo 
s Jo 



\Mt)\ 



(87) 



where rjs = \r]\vs, Kq is the modified Bessel function, and 
we have used the fact that 



d^q e'l"" 
"2^q2 + 1 



= ^o(|r|). 



(88) 



It depends on the encircling time T and the shape of the 
contour For a circular contour of radius R the inte- 
gral evaluates to TKq {R/£,)- For a large enough contour, 
then, ^dyn — exponentially over a distance ^. 
By contrast the topological phase 



7top = l^??l??2 



R"''^ It 



I X dx 



(89) 



does not depend on the encircling time T. For a cir- 
cular contour of radius R, the integral evaluates to 
27r[l - {R/^)Ki{R/0]- So, for large R 



7ri^?7i?72, 



(90) 



with an exponentially vanishing correction over a dis- 
tance ^. 



Ses ^ TTi^M' / d^fc j„^(-fc)M-;(fc)j:(^), (83) 



where 



1b = T^VV 



(84) 



By substituting the vortex current terms in Eq. (|83p we 
find the exchange Berry's phase 



hl^ / d^fc 5j'i{-k)M^^{k)53^2{k). (85) 



VI. EXPERIMENTAL SIGNATURES 

Fractional charge has been unambiguously detected in 
FQH systems using a quantum antidot electrometer— 
and shot noise analysis. ^^'^^ In other systems, such as 
the dimerized polyacetylene chain, ^ the experimental de- 
tection presents a greater challenge and the interpreta- 
tion of the results is less straightforward.— Experimental 
detection of the fractional statistics in FQH systems to 
date remains elusive. Although the claim has been made 
to detect the (abelian) fractional statistics in a Laughlin 
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stat#^ the interpretation of these results remains con- 
troversial. It appears difhcult to disentangle the effects 
of fractional charge from the statistics, these being, in 
essence, two complementary manifestation of the same 
underlying many-body wavefunction. It has been pointed 
out^i^ that the detection of the non-abelian exchange 
structure expected to occur in the Moore-Read pfafhan 
state could be in fact more straightforward, as this ef- 
fect is not directly tied to fractional charge. In general, 
unambiguaus detection of the fractional statistics in any 
physical system remains an unsolved problem and a chal- 
lenge. 

Against this backdrop we now discuss possible experi- 
mental signatures of the fractional charge and statistics in 
our proposed system. We outline a concrete experimen- 
tal setup for probing the fractional charge bound to the 
vortex defect and show how it can be measured directly 
in a bulk transport measurement. In fact this measure- 
ment appears to us more straightforward than any other 
scheme for fractional charge detection discussed in the 
existing literature. This simplicity arises from the fact 
that the number of vortices traversing the width of the 
system can be counted precisely through the Josephson 
relation, a fundamental property of the superconductor. 
The charge bound to these vortices, on the other hand, 
can be measured accurately owing to the precisely quan- 
tized Hall conductance, a fundamental property of the 
2DEG. 

Similar considerations allow for controlled manipula- 
tion of vortex defects. This opens up a possibility of 
moving them, in principle at least, along any desired tra- 
jectory in the system. We describe below an idea for the 
anyon shuttle, an all-electric system for anyon manipula- 
tion. The ability to move anyons in a controlled fashion 
should aid the future experiments aimed at probing their 
exchange phase, although drawing on the experience with 
FQH systems we expect this to present a significant ex- 
perimental challenge. 



A. Fractional charge 

The setup geometry is shown in Fig. [9] The system is 
subject to a perpendicular magnetic field |B| > -Bm, the 
matching field, which produces weakly pinned vortices in 
interstitial positions. In the 2DEG a charge 5Q = e/2 is 
bound to such defects. A supercurrent, JsCj is induced 
in the superconducting film producing a Magnus force, 
Fm = Jsc X B/c, on the vortices. The current density is 
largest in the constriction. By suitably tuning the mag- 
nitude of the current one could arrange for the force only 
to affect the interstitial vortices, generating a vortex cur- 
rent Ji, in the film along with an electric current Jg in 
the 2DEG as fractional charges bound to defects cross 
the constriction. The vortex and electric currents result 
in, respectively, a voltage drop, Vscj across the supercon- 
ducting film and a Hall voltage, V2deGj across the 2DEG. 
As we demonstrate below the ratio of these two voltages 




FIG. 9: (color online) Experimental setup for measuring the 
fractional charge. A supercurrent, Jsc, is induced in the su- 
perconducting film, which causes a vortex flow, J„, of un- 
pinned defects (white circles) on the background of pinned 
vortex lattice (black circles) across the constriction along with 
a corresponding electric current, Je, in the 2DEG. A Hall volt- 
age, V2DEG, is generated in the 2DEG and a voltage drop Vsc 
develops across the superconductor. This measures the frac- 
tional charge SQ — (V2DEG/2Vsc)e. The small arrows show 
the chiral edge current in the 2DEG and the thick (green) 
lines show the wires that carry the electric current in 2DEG. 



provides a direct measure of the fractional charge. 

The superconducting phase difference, At9, across the 
film has a time-dependence given by the Josephson rela- 
tion 



dAd 2e 



(91) 



dt h 

When a vortex crosses the constriction it accounts for a 
change of 27r in the phase Ai?, therefore 



1 dAd 2e 



2tt dt 



sc- 



(92) 



The corresponding electric current in the 2DEG is simply 
Je = SQ Jy. The Hall voltage is therefore found to be 



V2DEG — — |Je 



2SQ 



ve 



(93) 



where we have used the quantized value uh = ve^/h 
in the quantum Hall state. So, we arrive at a particu- 
larly simple relation for the fractional charge of defects, 
namely. 



5Q 



^2DEG . 



-ve. 



(94) 



The fractional charge bound to the vortex is thus sim- 
ply related to the ratio of two experimentally measur- 
able voltages. It is worth emphasizing that Eq. ([94]) is 
an exact relation whose validity relies only on the fun- 
damental properties of a superconductor and a 2DEG in 
the quantum Hall regime at filling fraction v. 
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FIG. 10: (color online) The "anyon shuttle." A superlattice 
of holes with diameter D and separation d, both ^ \l, is 
used as source and sink of anyons. The arrows on the lower 
left hole shows the chiral edge currents concealed on other 
holes for clarity. The black squares are the leads on the su- 
perconducting film. 



in various quantum computation schemes. An important 
feature of this scheme is its being completely electrical, 
foregoing the need for any moving parts. We believe this 
is a major advantage especially in view of possible appli- 
cations of non-Abelian anyons in quantum computation 
as briefly discussed in Sec. IVIII 

Using this shuttling scheme it should also be possible 
to perform bulk interference measurements to measure 
the fractional statistics of the anyons. A variety of such 
measurements exists, especially in the context of FQH 
physics. Assessing the merits of any of these setups in 
the current system, however, requires a detailed study 
and is beyond the scope of this work. We shall only 
note that an important condition for such interference 
experiments is the quantum coherence of the fractional 
bound states at least over the distance they move and 
the period of observation. While we do not see a general 
reason forbidding this coherence, attaining it in a real 
system will present a significant experimental challenge. 



B. "Anyon shuttle" and interference measurements 

Observing fractional statistics is a much more difficult 
feat requiring precise interference experiments. However, 
in our system, the fact that the anyons have a physical 
flux attached to them might prove useful. For instance, 
it is possible to perform hulk, as opposed to edge, mea- 
surements on the system. Bulk sources and drains of 
fractional particles may be created by building sources 
and sinks of flux tubes. Explicitly, we will in the follow- 
ing lay out an all electrical scheme for shuttling fractional 
particles around along any desired path of a superlattice 
of holes drilled in both the layers of the system. A sketch 
of one such design is shown in Fig. [TOl 

An array of holes drilled through both the supercon- 
ducting film and the 2DEG layer provides the sources 
and drains of anyons through a supply of superconduct- 
ing vortices trapped in the holes. The hole diameter, Z?, 
and spacing, d, between them are both much larger than 
the penetration depth A^. In equilibrium the magnetic 
fiux through each hole will be quantized in units of $o/2. 
In the 2DEG each hole has a set of low-energy chiral edge 
states with discrete spectrum due to the finite hole size. 
These edge states can accommodate extra charge which 
is crucial for the functioning of the anyon shuttle. 

Now imagine we attach wires to the superconducting 
film at the centers of the square leads as indicated in 
Fig llOl Through these wires we can feed appropriate su- 
percurrents. When the current is sufficiently strong the 
resulting Magnus force Fm moves the vortices. If the 
geometry is right vortices will move from one hole to an- 
other, carrying fractional charge with them. As before, 
the number of vortices traversing between the holes can 
be determined from the voltage drop via the Josephson 
relation, Eq. ([9^ . Using a suitable sequence of current 
pulses it is possible to shuttle vortices along any desired 
path in the device, effectively braiding them as required 



VII. CONCLUSION 

We showed in detail that the proximity of two well- 
understood, weakly interacting systems can result in 
some very interesting phenomena. Neither the super- 
conductor nor the integer quantum Hall state in the 
2DEG support fractionally charged quasiparticles. Their 
ground states can be understood by filling a set of single- 
particle states. However, when brought together, they 
give rise to excitations that carry fractional charge and 
obey fractional statistics. This result rests on the follow- 
ing general conditions: (i) the exact quantization of flux 
by the superconducting film in integer multiples of |<I>o, 
half the natural flux quantum of the 2DEG; (ii) the in- 
compressibility of the quantum Hall state of the 2DEG; 
and (iii) the pinning of the background Abrikosov vortex 
lattice by the artificial array of pinning sites. Conditions 
(i) and (ii) are manifestations of gauge invariance and are 
hence robust against weak perturbation (hence, disorder) 
in the system. Condition (iii) is necessary to keep the 
vortex lattice from adjusting itself when a quantized flux 
moves, basically allowing defects to be created. Other- 
wise, the incompressible vortex lattice will rearrange and 
thus screen any excess charge from accumulating in the 
2DEG. 

This general basis is what connects our different stud- 
ies in previous sections. The tight-binding model is rel- 
evant for strong Zeeman coupling due to large values of 
the gyromagnetic ratio." The continuum model of spin- 
polarized electrons with uniform background field is rel- 
evant for large values of penetration depth, Ib <C Al, 
whereas the intermediate regime was explored by the con- 
tinuum model of Pauli electrons with g — 2. Finally, the 
effective theory formulation allowed us to properly ac- 
count for the role of the superconductor and clarify the 
origin of the length scale ^, and in Sec. IIIII 

A conceptual subtlety must be noted in our usage of 
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the terms "quasiparticle" and "excitation" for the bound 
states produced by the vortex lattice defects. These 
bound states are clearly not excitations of the 2DEG 
in isolation. Indeed, inserting a flux ±^$0 changes the 
Hamiltonian, and hence the energy spectrum, of the 
2DEG. A true excitation of the isolated 2DEG is only 
found by inserting a flux ±$0, the natural flux quan- 
tum for unpaired electrons. However, these bound states 
are true excitations of the system as a whole in much 
the same way as the fractionally charged domain walls 
of the polyacetylene chain are true excitations of the ID 
electron-phonon systemi^ In both cases, the electronic 
system in the 2DEG or ID chain responds to defects in 
an ordered state arising from external interactions, in our 
case in the supercondcuting film and in the polyacetylene 
case in the ionic system. This is also clearly seen in the 
effective theory formulation of Sec. |Vl a vortex defect 
in the vortex-lattice phase is a finite-energy excitation of 
the whole system, shown to have fractional charge and 
statistics. 

We briefly comment on the range of parameters nec- 
essary to realize our proposed system in the laboratory. 
Several length scales must be considered: the penetra- 
tion depth Az,, the coherence length the magnetic 
length the thickness d of the superconducting film, 
the separation w between the superconductor and the 
2DEG, and the size L of the system. A natural order be- 
tween some of these scales must exist in order to realize 
the physics discussed in this paper. Specifically, we need 
£.0 d < Xl L. The requirement £,0 ^ ensures 
the superconductor is a strong type-II with high enough 
Hc2 ^ B the external field needed to realize the quantum 
Hall state in the 2DEG. In high-temperature cuprate su- 
perconductors Al ~ 1000 A and ^0 ~ 10 A and Hc2 can 
go beyond 100 T (Ref. S). Typical values of B are 5-15 
T in semiconductor heterostructures.— This corresponds 
to £b ^ 150 A, which determines the average spacing of 
the artificial lattice of pinning sites. This value is not 
too far from what is achievable in today's artificial struc- 
tures^Sii^ and can be increased if lower values of B are 
possible. Also, one could partially pin only a subset of 
the Abrikosov lattice using a larger lattice spacing with 
a specially engineered shape. For a large enough subset 
and the right shape, the rest of the vortex lattice will 
be pinned by vortex- vortex interactions. Defects may be 
created and moved the easiest when d ^ 0. However, the 
bulk penetration depth is modified in a thin film, given 
by the effective Pearl penetration depth Acff = A|^/(i, and 
diverges as d — > 0. This would destroy the bound state. 
Therefore the optimal conditions to allow the creation 
and motion of defects as well as the corresponding bound 
states are found when d < A^. We also require the size 
of the system to be larger than all other length scales to 
enhance the 2D nature of the system and minimize any 
field leakage from the edges of the system. Ideally, the 
separation w Xl so that most of the field that exists 
at the surface of the superconductor enters the 2DEG. 
However, we expect our conclusions to stand with minor 



corrections up to it» < Al. 

The free electrons with gyromagnetic ratio g = 2 con- 
sidered in Sec. lIIIBl mav not be just a theoretically conve- 
nient limit. Such electrons are realized on the surface of 
liquid helium and have been under investigation for some 
time.'*'' The typical electron density and mobility has so 
far been generally smaller than those in the semiconduc- 
tor heterostructures and, to the best of our knowledge, a 
quantum Hall state has not been obtained yet. 

Finally, we note that the same set-up can be used with 
the fractional quantum Hall states of the 2DEG. For in- 
stance, when V = 5/2, by the same general argument 
leading to Eq. ^ we find 



SQ5/2 



(95) 



This is the same charge carried by the true excitations 
of the Moore-Read state^^i^ thought to describe the 
ground state at 1/ = 5/2. This is consistent with the 
fact that the Moore-Read state can be thought of as a 
paired-state of electrons (corresponding to a p -I- ip super- 
conductor with a fixed particle number), which implies 
the natural flux quantum in this state is indeed 5$o- We 
then predict that the bound states found in such a set-up 
will behave as non-Abelian anyons.^" If so, our propos- 
als for detecting the fractional charge and the all-electric 
anyon shuttle in Sec. IVII gain new significance, as they 
can now be used to measure the fractional charge e/4 
and braid non-Abelian anyons. 
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APPENDIX A: SCHRODINGER EQUATION 
WITH (5-FUNCTION FLUX 

The Schrodinger equation TYpoi* = -E* for Eq. ^ is 
solved by a wavefunction ^ = e^^fR{r) in the dimen- 
sionless polar coordinate introduced in Sec. IIII Al In 
terms of the dimensionless energy e — 2E/hLu, where 



eBo/mgC is the cyclotron frequency (mei%uj — h), 



we find 



R" + -R' - 



-r"^ + 2(e + l-i-i) 



R = 0. (Al) 



Then by the usual Frobenius series expansion we can find 
the energy levels to be 



Jn) 



= 2fc+ 1^-77! - (^-77) + 1, 



(A2) 
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with the corresponding normalizable eigenstate 

fe 



(A3) 



and the coefficients 

s=0 



2(s - /c) 



(s + l)(2s + 2|/-r7| +3)' 



(A4) 



So for / > ?7 we find the levels E^j^'l = hw{k + i), 
which are the usual Landau levels. But for I < t] we have 
E^!^^ — huj{k ~ / + i) + rjhw, which are located inside the 
gaps between the Landau levels with a degeneracy k — 
I + 1. These states carry an angular momentum opposite 
to that of the usual Landau levels. This pushes their 
energy up but in the presence of the fractional extra flux 
this increase is not a whole integer multiple of the level 
spacing huj. We recover the full Landau level structure 
when 77 = or 1. 

The normalized states in the LLL are given by 



where z = x+iy is the (dimensionful) complex coordinate 
and z is the complex conjugate. This can be solved by 

V'o = f{z)e't (B4) 
where / is an entire function of z, d//dz = 0. 

For a uniform field, we have 



= (symmetric gauge), (B5) 

= T^iy^ '^^ (Landau gauge), (B6) 



Choosing f{z) = z'' in the symmetric gauge, we find at 
once the LLL states (77 = in Appendix [X]) . where I = 
0, 1, • • • indexes the angular momentum. 



APPENDIX C: TKNN INVARIANT FOR THE 
DIRAC HAMILTONIAN 



*[^)(z)- 



\z\ ''z'e 2 



v/27r£|r(l + ?-77)^ 



(A5) The negative-energy eigenstate of the massive Dirac 

Hamiltonian (1581) can be written as 



where the complex coordinate z — re^^ . 

APPENDIX B: AHARONOV-CASHER 
SOLUTION 

The Hamiltonian (l23ll can be written 



free 



2me 



(Bl) 



where 11 = p — |A is the dynamical momentum op- 
erator, and we have used the commutation relation 
[IIjjjIIj,] — \ehB/c. The ground states are zero energy 
states, TificcV'o = found by solving {ux^x + o-yRy)'<Po = 
0. Obviously ■00 rnust be an eigenstate of cr^. Renormaliz- 
able solutions only exist when the spin is aligned with the 
magnetic field, i.e. Uzipo = ipo, as physically expected. In 
the Coulomb gauge V • A = 0, we may write the vector 



potential in terms of a scalar 
which satisfies 

$0 



The zero-mode equation takes the form, 
dz az 



(B2) 



(B3) 



V2I ^l + m/E^ 



(CI) 



where 0p — (px — ipy)/\p\. Explicit computation using 
Eq. dniD gives 



1 z X p 
P " 2Ep{m + Ep)' 



(C2) 



It is easiest to evaluate Eq. (l56l) using the Stokes theorem, 



^ / d^p (Vp X ^p) 



(C3) 



27riBz'^ 2EI 



m 1 . N 
-sgn(m). 



where, in the last step, we have extended the upper 
bound of the integral to infinity. The latter approxima- 
tion is accurate as long as |m| <C W (the bandwidth), 
and becomes exact in the limit m — > when the Berry 
flux becomes a delta function at the node. Since /C is 
a topological invariant, constrained to be an integer, it 
cannot change as m is varied as long as the gap remains 
open. The above calculation thus gives the exact result. 
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